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It has recently been proposed in Refs. Q, Q that the dark energy could be attributed to 
the cosmological properties of a scalar field with a non-standard dispersion relation that decreases 
exponentially at wave-numbers larger than Planck scale (fcphys > Mpj). In this scenario, the energy 
density stored in the modes of trans-Planckian wave-numbers but sub-Hubble frequencies produced 
by amplification of the vacuum quantum fluctuations would account naturally for the dark energy. 
The present article examines this model in detail and shows step by step that it does not work. 
In particular, we show that this model cannot make definite predictions since there is no well- 
defined vacuum state in the region of wave-numbers considered, hence the initial data cannot be 
specified unambiguously. We also show that for most choices of initial data this scenario implies the 
production of a large amount of energy density (of order ~ Mp t ) for modes with momenta ~ Mpi, 
far in excess of the background energy density. We evaluate the amount of fine-tuning in the initial 
data necessary to avoid this back-reaction problem and find it is of order H/Mp\. We also argue 
that the equation of state of the trans-Planckian modes is not vacuum-like. Therefore this model 
does not provide a suitable explanation for the dark energy. 



I. INTRODUCTION 

Recently, it has been claimed in a series of articles 0, 0, 0] that the cosmic dark energy component could be 
explained naturally by the trans-Planckian energy of a scalar field with a suitable non-linear dispersion relation in the 
trans-Planckian regime. Such dispersion relations, which relate the frequency w p hy S to the wave-number fe p hys of a 
scalar field wave-packet, and which depart from the standard linear dispersion relation in the trans-Planckian regime 
are a way of modeling phenomenologically the unknown physics for sub-Planckian wavelengths. They have been used 
extensively in the recent literature in the context of black-hole physics Q and of the inflationary trans-Planckian 
problem |5|. 

In the particular case considered in Refs. 0,0,0, the dispersion relation departs from its standard linear form 
and approach a decreasing exponential at large wave-numbers. This type of dispersion relation could possibly emerge 
from string theory . It has been argued that the energy density of the modes of sub-Planckian wavelengths and 
sub-Hubble frequencies (referred to as "tail" modes) is naturally of the same order as the critical energy density today 
and has the same equation of state as a cosmological constant. Hence, it could account without fine-tuning for the 
dark energy. The energy density contained in the tail today has been calculated in Ref. by solving for the time 
evolution of a test quantum scalar field evolving in the curved cosmological background, assuming its initial state 
at the onset of inflation is the vacuum. The equation of state of the tail modes has been calculated in Ref. [3| and 
its cosmological evolution has been solved to argue that the cosmic coincidence problem (why does the dark energy 
dominate now?) is solved. 

In this paper we argue that this model does not and can not work for several reasons. We first argue that the energy 
density and equation of state of the "tail" modes depend directly on the choice of initial data for the scalar field, 
and that this latter cannot be specified unambiguously since there is no preferred initial vacuum state (Section [Hj . 
This implies that any cosmological consequence derived depends directly on the ad-hoc choice of initial data (initial 
quantum state). We then show that the violation of the Wentzel-Kramers-Brillouin (WKB) approximation in the 
remote past for all co-moving wave-numbers, which is inevitable in the present scenario, implies the continuous 
production at all times of a large amount of quanta with physical wave-number ~ Mp\ (Section . This finding is 
in agreement with general arguments given by Starobinsky Q (this latter work did not however study the present 
scenario). We evaluate the amount of energy density produced for modes of physical wave-number and frequency 
~ Mpi for various choices of initial data and conclude that it is generically of order Afpj. This process takes place 
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FIG. 1: Dispersion relation u> p hys vs fc p h y s- Region I (hashed area) corresponds to the "tail" modes for which < H [and 

fcphys > K(rfj\. Region II corresponds to the sub-horizon modes (fcphys > H) that are outside of the tail, while region III 
corresponds to the super-horizon modes k p h ys < H. The fundamental scale k c is also indicated. 



at all times, and since the energy density produced is much larger than the background energy density, it implies 
that the semi-classical perturbative framework on which the model of Refs. 0, 0, 0] rests breaks down. In an earlier 
study, devoted to constructing an effective stress-energy tensor for theories with non- linear dispersion relations ja], 
we already criticized this model by arguing that it led generically to the wrong equation of state. We revisit this issue 
further in Section llVI where we prove that the effective energy-momentum tensor we derived earlier is well-behaved, 
thus disproving an improper claim of Ref. Q and confirming our earlier criticisms. We provide a summary of our 
conclusions in Section Ivl 



II. INITIAL CONDITIONS FOR THE MODE EVOLUTION 

Refs. ! i consider a scalar field, <fi, with a non-linear dispersion relation that is linear in the sub-Planckian 
regime and approaches a decreasing exponential at trans-Planckian wave- numbers (i.e. for fc p h ys ^ k c , k c ~ Mp\ 
being a fundamental characteristic scale). This dispersion relation is shown in Fig. ^ This scalar field is assumed 
to describe the density (scalar) perturbations and/or the primordial gravitational waves. The "tail" modes are thus 
interpreted as a bath of gravitons of super-Planckian wavelengths and sub-Hubble frequencies. This scalar field is 
treated as a test-field (its back-reaction on the background is neglected) and is quantized on the curved cosmological 
background. Assuming that the "tail" modes of this field are initially in a well chosen vacuum state as 77 — > —00 (rj 
denoting conformal time), the occupation number at late times (rj — > +00) of quanta extracted out of the vacuum 
by the dynamical background has been calculated in Ref. 1] . This occupation number can then be used to calculate 
the energy density stored today in the "tail" . This is the thread of the calculation performed in Ref. [Ij , which we 
now follow in some detail. This discussion will take us to the two main arguments that we bring forward against this 
model (given in this Section and the following). 

The equation of motion of a scalar field = /i/a in a Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time 
with scale factor a{rf) reads: 



/4' 



^-(1-60- 
a 



Mfc = 0, (1) 



where £ is a coupling parameter to gravity and a prime denotes differentiation with respect to conformal time. 
k = afcphys is the co-moving wave-number and ui = acjphys is the co-moving frequency. £ = for tensor perturbations 
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degrees of freedom and £ = 1/6 for a conformally coupled field. The scale factor is taken to be a power law in 
conformal time, a = |fy/?7c| , and the following dispersion relation, parameterized by two parameters t\ and €3 with 
£3 = 4 — 2ei in order to insure that the dispersion relation is linear for small wave-numbers, is introduced: 



w phys(^phys) — ^phys 



6i 636" 



(1 + e*) 2 



(2) 



with x = (fcphys/^c) 1 ^ = A\r}\ with A = (k/kc) 1 ^ /\r] c \. A problem with this dispersion relation is that it depends on 
the power-law index [3 of the scale factor. If taken literally, this means that the dispersion relation, or the physical 
frequency of a given mode, changes as the scale factor power-law index (3 changes between various cosmological 
eras (e.g. inflation / radiation domination / matter domination). More importantly one easily sees that the above 
dispersion relation has a pathological behavior in the radiation [(3 = — 1) or matter ((3 — —2) dominated eras. In fact, 
for (3 < 0, it implies Wphys/^phys — > as fc p h ys — ► 0, whereas one should instead reach the linear dispersion relation 
in that regime with ui p hys/k p h ys — > 1. Since Ref. [j focused on the case of de Sitter space-time with (3 = 1, we set 
f3 = 1 in the above dispersion relation, i.e., the above x should be understood as x = k p hy S /k c . This reformulated 
dispersion relation thus coincides with that used in Ref. Q for de Sitter space. However, in the matter dominated 
era, for instance, we have x = fc p hys/fcc °c r\~ 2 and the general class of solutions to the field equation obtained in 
Ref. does not hold anymore. The linear dependence of x on 77 is lost for background metrics other than de Sitter, 
but the linear dependence of w p h ys on A; p h ys is preserved in the small fc p h ys limit for all metrics, which is obviously an 
imperative. The field equation can finally be rewritten as: 



ei e 3 e" 



(l + e*) 2 



with, again, x(rj) = fc p h ys /fcc = k/[a(r])k c }. The solution to this equation depends on the value of £ and (3. In Ref. Q, 
the contribution of the [a" /a] term is assumed to be negligible at early times. However the above equation shows that 
this is not the case; denoting by ft 2 the term in curly brackets in Eq. ©, one has fl 2 ~ k 2 (ei + e 3 )e~ k ^ f / C fe <=l r 7<=l /3 ) — 
(3{f3 + 1)(1 — 6£)/?7 2 as 77 — » —00 and the term a" /a oc vT 2 is always dominant in that limit if £ 7^ 1/6 ((3 > 1, 
77 — > —00). Therefore, in the limit £ 7^ 1/6 and 77 — > —00 the two independent solutions to the field equation are 
power-laws in 77: 



> «± = ^±^+/3(l + /3)(l-6£). (4) 

This is an important point as it implies that the mode function does not behave as a plane wave in the limit 7/^—00 
when £ 1/6. The solution to the field equation in this limit is reminiscent of the mode freezing in inflationary 
theories for fields with linear dispersion relations and £ = when the mode exits the horizon. 

It is also argued that the term in a" /a can be absorbed at late times in a redefinition of the dispersion relation. 
However this cannot be correct since by construction, the dispersion relation w p h ys (fcphys) depends only on fc p h yS j i-e. 
time only enters via fc p h ys - Therefore, one can absorb a term a" /a oc r/~ 2 into ^ p h ys / a2 onr y if /3 = 1 (de Sitter 
space-time), as inspection of Eq. Q reveals. In effect, the curly bracket of Eq. J2Jl can then be rewritten as or 2 
times a function of kr\ ex fc p hys- But, in that particular case, the redefined modified dispersion relation does not have 
an exponential shape anymore, since u> p h ys approaches a constant (~ H) as fc p h ys — ► +00. However, this should not 
give the impression that the corresponding solution is a plane wave since, evidently, the co-moving frequency u> which 
enters Eq. Q still behaves as oc r\~ 2 . Moreover in the case of a matter or radiation dominated cosmology, one cannot 
absorb the scale factor term in the dispersion relation. 

Nevertheless one can also assume £ = 1/6. In that case, it is possible to find an exact solution to the equation 
of motion in de Sitter space-time. Indeed, for a conformally coupled field, the term a" /a disappears from the 
field equation and the equation becomes simpler. Note however that the scalar field cannot correspond to tensor 
perturbations degrees of freedom since these arc minimally coupled to the metric. Let us consider £ = 1/6 for the 
moment. A solution to the field equation, given in Ref. pj reads 1 : 



/4 m) (>7) = C (in) (1 + e-*)' 
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d+ r- b 
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(5) 



1 Equation (25) in Ref. y| contains a misprint that has been corrected in the following equation 
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where 2 ^i is an hypergeometric function and b and d are expressed in terms of ei and £3 as: 



h, d =\/\ + i 3> and e 4 = ^c ? 7c e J- (6) 

This solution is valid only for de Sitter space-time with x = kr]/k c r] c (77 < 0, r) c < 0). As already mentioned, this is due 
to the fact that, with the reformulated dispersion introduced above, the linear dependence of x in the conformal time 
is lost for other scale factors. However similar solutions for other metrics can be obtained if the dispersion relation 
is tuned to the power-law evolution of the scale factor, i.e. if the parameter x remains linear in 77 (possibly at the 
expense of linearity of w p h ys in the small fc p h ys limit, see above). Equation Q has in fact two independent solutions 
(see below) and the choice ||5J represents only one branch of the solution, which is moreover written on the branch 
cut of the hypergeometric function 2-Fi- At early times (77 — > —00, i.e. x — > +00), this solution (JHJ does not oscillate 
and it blows up 2 . 

It is more convenient to write the general solution to the field equation, with £ = 1/6, as = Ci/xjj, + C^/J^, 
where jj,^ and /i? are two independent solutions given by: 

4%) = e bx (l+e x ) d + 1 /\F 1 (b+d+±, b+d+ ~; 2b +l;-e x ^j 

$\r,) = e- bx {l+e x ) d+1 l 2 2 F 1 ^-b + d+^-b + d+^-2b+l-,-e^j. (8) 

Since b is pure imaginary, and d is real, one concludes easily that /A = lA. ■ The Wronskian of these two solutions is 
non-zero, and can be used to relate the coefficients G\ and C 2 so as to obtain canonical commutation relations for the 
field operator and its adjoint. Since only one branch of the solution was given in Ref. the canonical commutation 
relations for the field and its adjoint could not be satisfied. More precisely, it can be checked that the solution given 
in Eq. j£| is real. This is due to the fact that it involves an hypergeometric function of the form 2-F1 (a, a*;a + a*; z) 
with z = z* — 1 + e~ x and a = b + d + 1/2 in that case and 2 F* (a, b; c; 2) = 2F* (b, a; c\ z) = 2 Fx(a* , b* ; c* ; z*) . 
Therefore, one has 2-Ff (a, a*; a + a*; 1 + e~ x ) = 2 Fi(a, a*;a + a*; 1 + e~ x ) and the mode function is indeed real. 
It follows that the Wronskian of the solution considered in Ref. Q vanishes: W(fi, fJ,*) = fikfJ-k' ~ Vh^k — 0- Using 
the properties of hypergeometric functions, one can check that both independent solutions fjjp and ju,^ behave as 
oc x oc \r\\ in the limit r\ — > —00, i.e. these mode functions blow up. This result is consistent with Eq. (JJJ) since, in the 

tail, the two branches and /u,^ are linear combinations of the Bessel functions Jo and Nq. 

Therefore, we have shown that neither in the case £ 7^ 1/6 nor in the case £ = 1/6 does the mode function behave 
as a plane wave in the tail. Thus the initial state of the field cannot reduce to the Bunch-Davies adiabatic vacuum, 
contrary to the claim Q: "we show that there is no ambiguity in the correct choice of the initial vacuum state. The 
only initial vacuum is the adiabatic vacuum obtained by the solution to the mode equation" . The usual prescription 
to remove the ambiguity on the choice of vacuum state in curved space-time, i.e. for constructing a vacuum state 
which is closest to the definition of vacuum in Minkowski, is indeed to rely on the WKB approximation to construct 
vacua of successively higher adiabatic order 0. In this scenario 0013, this construction cannot be performed 
for a simple reason: the WKB approximation, which quantifies the adiabaticity of the quantum mode evolution is 
violated at all times for modes contained in the tail, i.e. modes with fc p h ys > &c and w p h ys < H. The WKB condition 
can be written in the form |Q/ri 2 | <C 1 0]> where fi 2 , denotes the term in curly brackets in Eq. as before, and 
Q = f2'//2Slfc — (3/4)Sl^? /^|- In effect the WKB solution /xwkb = exp(±i J ttkdT)/^2Qk exactly verifies the following 
equation /^wkb + (^fc ~ Q)a*wkb = 0. Therefore it is a good approximation to the solution of the actual mode equation 
fi" + fi 2 /i = if the above inequality |Q/r2 2 | <g; 1 is satisfied (see also Ref. ^lj for more details). 



2 A hypergeometric function of the form %Fi (a, (3;a + /3; z) is singular at z = 1. 
limit r\ — > —00. In this case, the equation reduces to + fc 2 (ei + es)e~ x fi^ = 



One can also solve Eq. J3J for £ = 1/6 and /3 = 1 in the 
and the solution can be written as 



1 



Ai(fc)Jo(ef) + A 2 (k)N (ey) 



(7) 



where Jo and No are Bessel functions and where y = — x'/2 + ln[4fc 2 (ei +e$) /A 2 ]/2. The Neumann function diverges in the limit r\ ■ 
(y — » —00). In the tail, the corresponding behavior for the scalar field itself is given by <j> oc r\ and (j> oc rp . 
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The expression for Q/f2| is cumbersome, but since we are interested in the regime fc p hy S 3> k c , we may use the 
limiting form of the dispersion relation: 



w phys ~ fc phys Vei + £se fc pW( 2 M (fc phys > k c ), (9) 

If £ = 1/6, then = aw p hy S and |Q/fJ|| ~ k p \ lys .ff 2 /16feca; p ! hys . l n or der to understand the behavior of |Q/n||, it is 

convenient to introduce the physical wave-number K + > fc c such that uj p Y iys (K + ) = + 0)/0H [in the case £ = 0, 
one has ttk{K + ) = 0]. This wavenumber K + ss K, where K (77) is the physical wavenumber that gives the lower limit 
of the tail, as indicated in Fig.^ Using Eq. ©, one easily derives: 



K + ~ 2fc c ln 




(10) 



This formula is written to zeroth order in ln(fc c / H)/(k c / H) but can be expanded to arbitrary order in a straightforward 
way. The meaning of the physical wave- number K+ is the following (see Fig. If fc p h ys <C K+ but fcphys > fcc 
(i.e., within region II), the mode is outside of the tail with w p h ys 3> H. If however fc p h ys 3> K+, the mode is in 
the tail with uj p hy S -C H (region I of Fig. P). Then fc p h ys S> K + means that co 2 hys <C H 2 which implies in turn 
|Q/0^| ^> fcph ys /(16fcc) ^S> 1, hence the WKB approximation is violated at all times in the tail. Note that outside of 
the tail, i.e. for fc p h ys <C K + and fc p h ys 3> fc c (region II of Fig.QJ, the WKB approximation becomes valid. In effect 
fcphys/^phys a exp(fc phys /2fc c ) < \n(k c /H), hence |Q/fi|| < 1. 

If £ ^ 1/6, then for w p h ys if (region I or tail in Fig. QJ, the dominant term is cx a" /a in the expression of CI 2 ., 
namely - -(1 - 6£)(1 + 0)a 2 H 2 /0, hence - [4(1 - 6^)0(0 + l)]" 1 , which for = 1 (de Sitter) and f = 

(minimal coupling) reduces to 1/8. In this case, it can be shown that the WKB^ approximation does not give the right 
behavior for the mode function even though |Q/0|| is smaller than unity ^l|]> and that the WKB approximation 
is not valid either. Again, note that outside of the tail (region II of Fig. Q the WKB approximation is valid. The 
calculation is the same as in the previous paragraph, since for fc p h ys <C K + and fc p h ys fcc, one has fl 2 ~ a 2 to 2 hys 
since w p h ys 3> H. Thus one hnds |Q/n?| <C 1 outside the tail even for ^ 7^ 1/6. 

To summarize this discussion the WKB condition is violated by the present dispersion relation in the tail (region I 
in Fig. ^) at all times and an initial vacuum state cannot be constructed unambiguously. Outside of the tail (region 
II of Fig. Iy, for ojphys 3> H or fc p h ys *C -K+, the WKB approximation is a good approximation. One can also verify 
that the construction of an initial vacuum state by minimization of the energycontent does not work in this case, 
see Ref. 0. This point is one major obstacle to the scenario proposed in Ref. Since there is no preferred initial 
vacuum state, all cosmological conclusions drawn depend directly on the particular choice of the initial state, hence 
on the choice of initial data. At the very best, one has to fine-tune the initial conditions to obtain a given amount of 
energy in a given part of the spectrum. 

The standard calculation of the amount of energy contained at late times in a given co-moving wave-number mode 
is done by decomposing the solution at late times (outside the tail) in terms of positive and negative frequency plane 
waves, as 



07, 
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(ii) 



The squared modulus of the Bogoliubov coefficient 0k then will give the occupation number of quanta produced in the 
mode of co-moving wave-number fc. Note that, in principle, the coefficients ctk and 0k can be slowly varying functions 
of time, and the above expression implicitly involves a WKB approximation to first order in which the time evolution 
of ckfe and 0u is neglected. The corresponding vacuum is an adiabatic vacuum to first order. 

In Ref. [l| 0k is calculated in the limit r\ — > +00 as oj out — ► ^/eik. However the limit rj — > +00 does not hold in an 
inflationary Universe with a cx |?7|~' 3 and > 1 since a is singular as 77 — ► 0~. One needs to match the background 
evolution to a decelerated Universe as ry — > CP. In effect, if one wishes to calculate the contribution of the tail modes 
to the energy density today, it is necessary to calculate the evolution of the modes from the inflationary era up to 
today. Note that the dynamical evolution of the tail modes depends a priori strongly on the background scale factor 
dynamics. 

This calculation could not be performed in Ref. [l| , since the solution given in terms of the hypergeometric function 
is not valid at late times in the radiation dominated or matter dominated eras unless the parameter x of the dispersion 
relation is tuned to the evolution of the scale factor, but the dispersion relation would become pathological as we saw 



6 



before for x = (fephys/^c) 1 ^ with (3 < 0. Furthermore, as explained above, the solution to the field equation given 
in Ref. [see Eq. describes only one branch of the solution. Finally, one cannot compute /?/. for modes still 
contained in the "tail" at late times by matching the solution to plane waves as done in Ref. [Ij since for those modes, 
the WKB approximation is never valid so that the out solution cannot be decomposed in a sum of plane waves. Thus 
the calculation of the Bogoliubov coefficient (3k performed in Ref. [l| cannot apply to modes contained in the "tail" 
today. 



III. THE TAIL ENERGY DENSITY 



In this Section we calculate the amount of energy density created in quanta that redshift out of the "tail" , and 
show that it leads to a severe back-reaction problem. In Ref. pjthe energy density contained in the tail is calculated 
as 

1 f +co f 

(Ptail) = fcphysdfcphys J UJ phys duJ phys \(3 kphys | 2 , (12) 

where kn is the physical wave-number such that wh = Wphys(fcphys) = Ho today. This expression for (ptail) is ill— 
defined due to the double integration element dfcdw in the absence of a Dirac function on the mass shell. The 
total energy density (ptotai) is defined analogously but the lower bound is extended to k — 0. Then, it is argued 
that (ptaii)/(ptotai) — 10~ 122 during inflation. Note that if p ta n — 10~ 122 Mp[ and (ptaii) is constant (corresponding 
to a acuum-like equation of state as suggested) in order to account for the dark energy then the above statement 
yields ptotai ~ -^pi- If this holds during inflation, one faces a severe back-reaction problem since the background 
energy density during inflation is ~ 10 orders of magnitude below Mp,, and the overall calculation framework (a test 
quantum scalar field on a classical background) breaks down. As we argue in this Section, it is actually a generic 
prediction of this model that ptotai ~ Mpi a t a ^ times. This result ptotai ~ -Mp] is in agreement with a recent work by 
Starobinsky which showed that models with dispersion relations such that the WKB approximation is not valid 
in the far past when the physical wave-number fc p h ys 3> -Wpj implies a very substantial amount of particle production. 

In the following we calculate the amount of energy density stored in modes with physical wave- number fc p h ys ~ M-p\. 
The calculation follows the line of thought indicated in the previous Section. Since in the range H <C fc p hys <C K+ 
the WKB approximation is valid, one can decompose the solution to the field equation in terms of plane waves as in 
Eq. I|ll|) when modes enter this regime. As long as |Q/fi|| <§C 1 one can neglect the time evolution of (3k, and it is 
natural to interpret \(3k\ 2 as the occupation number of particles in mode k. As argued earlier this decomposition in 
plane waves cannot be made for modes that are still contained in the tail. 

One can then calculate the amount of energy density dp w /dln(fc p h ys ) stored in the log interval around the physical 
wave-number fc p h ys and the corresponding fractional density parameter df2 w /dln(fc p h ys ) in units of the background 
energy density: 



d ^ = 4 kphys^phys 1 2 ,,o\ 

dlnfc phys 3^ IP A/*, IPfe| ' 1 J 

using dp^/d ln(fc p hys) = fc p hys w phys|/^fc| 2 /2' 7r2 - The fractional density parameter must be smaller than unity at all times 
and for all physical wave-numbers, otherwise back-reaction is significant and all semi-classical first order calculations 
arc unreliable. In the following we calculate this quantity df2 w /dln(fc p hy S ) for a physical wave-number fc p h ys ~ k c , i.e., 
once the wavelength becomes larger than the fundamental scale. It can be expressed via (3k in terms of the constants 
that parametrize the choice of initial data. Our goal here is to study the dependence of the amount of energy density 
created in modes of physical momenta ~ Mpi on the initial data, for which there is no definite prescription as we 
argued in the previous Section. 



A. Conformal coupling: £ = 1/6 

In the case of conformal coupling £ = 1/6 there exists an exact solution to the field equation written in terms of 
the two independent solutions fJ^' and ^ in Eq. ©. One can then calculate the Bogoliubov coefficient deep in the 
region where the WKB approximation is valid, for instance around fc p hy S ~ &c by decomposing this exact solution in 
plane waves. However the coefficients of the hypergeometric function in term of which the exact solution hence (3k 
are written are of order k c /H 3> 1. For values of these coefficients that are relevant for our cosmological applications 
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(i.e. k c /H ~ 10 6 during inflation), the numerical calculation of the hypergeometric function turns out to be too 
involved and we have been unable to calculate (3k in a reasonable amount of time for k c /H > 10 3 . Therefore we take 
a different approach and approximate the exact solution in the tail fc p h ys > K+ by the solution derived in terms of 
Bessel functions in Eq. Q, and that in the region k c < /c p h ys < K + by the plane wave solution. The Bogoliubov 
coefficient (3k of the plane wave solution is obtained by matching the two solutions and their first derivatives at the 
transition point fc p h y s = K+. Of course, it gives an approximation to (3k, but as we show in the following the deviation 
from the overall behavior of (3k away from its minimum and on its behavior around its minimum are negligible. We 
thus proceed as follows: in the following Section, we calculate the Bogoliubov coefficient denoted ^j- approx -' by solving 
for the Bessel functions in the remote past and performing the matching at K + . In the subsequent Section, we 
calculate the Bogoliubov coefficient ^( exact ) analytically using the exact solution and demonstrate that j g( approx ^ i s a 

good approximation for values of k c /H as high as ~ 10 3 . Finally we examine the behavior of j g( approx ) and evaluate 
the amount of energy density produced by the non-adiabatic evolution of modes in the "tail" for realistic values of 
k c /H. This calculation is entirely analytical, only the verification of the accuracy of the approximation is numerical. 



1. Approximate calculation of the Bogoliubov coefficient 



As already mentioned above, see Eq. J7J, the mode function in the tail can be approximatively expressed in terms 
of the Bessel and Neumann functions Jo and Nq as 



Mfc(»7) 



1 



Ai(fc)J (e») + A 2 (k)N ( e y) 



H r !, 

v = kri H — in 

' J 2k r 2 



4fe c 2 (4- £l ) 
H 2 



(14) 



This solution is valid if the scale factor is that of the de Sitter space-time: a(ij) = —l/(Hr)). The mode function 
must satisfy the relation W = MfcMfc' ~ MfeMfc — *■ Using the above equation, one finds that the Wronskian is equal 
to W = —H(A2A\ — A\Al£) / {2Ttk c ) . As a consequence, if one represents the coefficient A2 in polar form, A2 = re 1 *, 
one has A\ = — nk c /(Hr sin $), where we have chosen A\ to be real. The parameters r and $ will characterize in the 
following the choice of initial data. 

In the region where the WKB approximation is valid, i.e. for w p h ys ^ H, one has 



Oik 



„-in 



VMM) 



VMk~n) 



(15) 



where fl = f 1 d.TUj(k,T). In order to express the Bogoliubov coefficient \(3h\ in terms of the constants parameterizing 
the choice of the initial data in the tail, A\{k) and A^^k), we use the continuity of the mode function fik and of its 
derivative at the transition between the two regions at y — y nl , for which o; p hys(ym) = \[2H . The result reads 



o(approx) _ 

Pk — 



IC 



7 fc (2/m)Jo(e^) + f |"e«°Ji(e»" 



y/4kuj(k,y m ) 

where 7^ = u>'/(2u>) + iu. Working out this last expression, one obtains 



+ A 2 (k) 



■7*JVo(e 1 "») + ~e»»JVi(e«») 



(16) 



4' 



(approx) 



( H 



yfW2\ K ^ 



1/2 



T~ <Mk) 



- ( ^+iV2 ) Jo(e^)+^i(e y ") 



+A 2 (k) 



V2 ) A (e^ )+V2N 1 (e y - ) 
4fc c J 

(17) 

We are now in a position where we can compute the | / g^ approx ) | 2 using the parametrization of the coefficients A\ and 
A9 introduced above. The final result reads 



|/3£ approx) | 2 (p,$) 



7T 2 J Np 2 nK 1 
1=— 7, o 1 7= =COt$ — - 

4V2p 2 sin 2 $ 4V2 2V2 2 



(18) 



where we have defined the rescaled variable p by p = ry H/K+ and where the coefficients J, A and K can be 
expressed as 



J = 
K = 



1G 

K. 



J__fcc 

V2A 4 



J Ji+2^(J 2 + J 2 ), N 



A' 2 



1 K. 



16/c c 2 V2 K 



■A Ai + 2(iV 2 + Nl) , 



^oAo + ^(JoAr + JrTVo) + 2-^(J N + JiAr) 
1ok c 4 K + 



(19) 
(20) 
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The Bessel and Neumann functions are evaluated at the matching point for which their argument reads e ym = 
2\/2kc/K + . A direct calculation shows that JN — K 2 = 2/tt 2 . The Bogoliubov coefficient | / #( approx J | 2 can b e viewed 
as a two-dimensional surface parametrized by the polar coordinates (p, <£>). 



2. Test of the method of approximation 

Before studying the above Bogoliubov coefficient in greater detail, one must check that the approximation is well- 
controlled. For this purpose, it is interesting to calculate the Bogoliubov coefficient using the exact solution expressed 
in terms of hypergeometric functions 



1 



(21) 



where the functions fi^ and /i^ have been defined in Eq. (JSJ) above, and the (dimensionless) functions Ci(fc) and 
C 2 (k) are related to each other by the Wronskian normalization condition [12^ : 



w4'-/4/4 = [|Cx(fc)| 2 - |c 2 (fc)p 



2fc c |?7o 

with the numerical factor T is written in terms of the parameters b and d, as: 



T. 



(22) 



■d+1) 



2b +1 



2b 2 F 1 
(-b + 



,F\ ( b + d+~,b-d+~;2b- 



1 , , f 



2;-l 2*1 



-b+d- 



b + d+ -, b - d + -; 2b + 1; -1 2 *i [ -b - d + 



2' 



-b + d- 



2b- 1 



F x [b 



d+-,b-d+-;2b 
2 2' 



2*1 



-;-26+l;-l 



-26+1; -1 



-b-d + 



-b + d 



-26 + 2:-l 



(23) 



This solution is valid at all times since it is an exact solution of the field equation. In this case, one can calculate the 
Bogoliubov coefficient at any time provided the WKB approximation is then valid, using: 



(exact) 



l 



2^ 



(24) 



where, in the last expression, is given by Eq. Q21[l. Notice that this procedure differs from the previous calculation 
of the Bogoliubov coefficient. Here, we do not perform a matching at the transition between the tail and the WKB 
region but rather use the exact solution (|2~T)> all the way through and calculate its "overlap" with the WKB solution 
deep in the WKB region. The initial conditions enter this expression via the two constants Ci(fc) and Ci(k). 

We need to compare |/?^, cxact ^| with |/3^. apprtB: ^| for the same initial conditions. Since |/3^ appraK ^| is expressed in terms 
of the constants Ai(k) and ^(fc), one needs to re-express A\{k) and A 2 {k) in terms of C\{k) and C2(k). This can 
be done by matching the asymptotic behaviors of the two solutions deep in the tail, i.e. in the limit r\ — > — oo. There, 
the approximate solution given by Eq. (|f 4f> reduces to 



2k\ 



nk c 



2A 2 , 2A 2 A 2 , 
Iln2 + l 7E + ^ln 

7T 7T 7T 



4fc c 2 (4 - ei) 



H 2 



(25) 



where 7e is the Euler constant, 7e — 0.5772. On the other hand, the exact solution of Eq. (|2f I) can be written as 



1 H 

2k k c 



k\ V \\ (GCi + G*C 2 ) - Gd 



27e + * [b + d 



V\b-d 



-G*C 2 



27e + * \-b + d 



Vl-b-d 



(26) 



where the coefficient G is given in terms of the Euler Beta function as: G = \/B(b + d + 1/2, b — d + 1/2) = 
r(26+ l)/[T(b + d+ l/2)r(&- d+ 1/2)], and satisfies, since b is pure imaginary and d is real, G* = G(b <-> -6). This 
relation stems from the asymptotic behavior of the hypergeometric functions for large values of their argument given 



FIG. 2: Left panel: The function A defined in Eq. 12911 that quantifies the deviation of the approximate \/3k\ with respect to its 
exact form. A is depicted as a function of (R,8) that parameterize the initial conditions (Ci,C2) of the exact solution. Right 
panel: zoom around the origin in order to unveil the two peaks structure. 



by Eq. (15.3.13) of Ref. Qj. The di-Gamma function ^(x) function is defined by ^(x) = dlnr(x)/da;. If one identifies 
the constant term and the linear term in conformal time of the two previous relations, we obtain: 



A 1 (k) = Gd(k)hn 



rfc c 2 (4- £l ) 



+G*C 2 {k)\]n 



H 2 

k 2 (4 



* b + d 



EP 

A 2 (k) = -n[GC 1 (k) + G*C 2 (k)} 



-b + d + 



* b-d 



1 
2 

-b-d + 



(27) 
(28) 



Then, it is sufficient to use the above relations in Eq. I|17l) to obtain \(3 t 



(approx) I 



in terms of C\(k) and C2{k) and 



compare it to |/3^ cxact ^|. In order to characterize the accuracy with which the Bogoliubov coefficient is calculated, we 
plot the following quantity 



A = 2 



|^(approx) | 



(exact) 



|^(approx)| |^(cxact) 



(29) 



for various values of the coefficients Ci(fc) and C 2 (k). More precisely, we use a polar representation and take C 2 (k) = 
Re* 8 while Ci(fc) is real and calculated in terms of C 2 (k) using the Wronskian relation Eq. (|22|l . In Fig. [3 we have 
plotted A(R, 9) for k c /H = 10 2 and k c = Mpi. We see that the error for large values of p is less than ~ 40 % and is 
constant, i.e. the offset between the two Bogoliubov coefficients does not depend on p and $ in a first approximation. 
For p ~ 0, the error increases to 1; this artefact comes from the fact that the minima where the two Bogoliubov 
coefficients vanish are slightly offset one from the other. If one coefficient vanishes while the other remains finite and 
non-zero, then the value of A is pushed toward one, A — *■ 1. This error however is of no consequence for what follows. 
Indeed we will not be interested in the location of the minimum but in the behavior of (3k around the minimum and 
far from the minimum. As is obvious from Fig. El these behaviors match closely in both cases and our approximation 
will be sufficient for our purposes. We have checked that the function A remains the same for other values of k c /H 
which allow numerical calculations, i.e. k c /H £ [10, 10 3 ]. 

The situation is in fact very similar to the standard calculation of the power spectrum in an inflationary theory: 
in principle, one cannot match two different branches at a point where the approximation breaks down (for standard 
inflation this occurs at first horizon crossing). However, since the approximation is only violated in a small region 
one expects the corresponding result to be correct at leading order. This is indeed the case for inflation for which 
the amplitude of the spectrum is predicted up to a factor of order unity and the spectral slope is unchanged by the 
matching. Here we also find that |/?[ approx) | = C(l)|^ cxact) |. 
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3. Fine-tuning of the initial conditions 

Since we have demonstrated that the approximation to \(3k\ is quite reasonable, we now study the behavior of 
Eq. I|18|) for more realistic values of the ratio k c /H. The Bogoliubov coefficient possesses an absolute minimum with 
/3fe = located at 



4 

Pmin 



7T J 



JN — K 2 



cos <£> r 



K 



(30) 



using the notations defined previously. One should not be surprised to find a minimum with 0k = since one can 
express the matching of the two branches of the solution and their first derivatives at rj m as two equations relating 
the coefficients A\(k) and A^ik) as a function of a>k and (3k and find a solution with (3k = 0. The Wronskian 
normalization condition is always satisfied by both branches of the solution. This solution with initial conditions 
(Pmin, 3?min) corresponds to a choice of initial data such that at late times, when modes have exited from the tail, 
their quantum state is that of an adiabatic vacuum. Note therefore that there is no naturalness in choosing these 
initial conditions since the adiabatic vacuum is only a late time consequence of such initial data. Furthermore one 
can show that for generic initial data, the state of the quantum field at late times is not an adiabatic vacuum, hence 
quanta have been produced. 

Indeed the behavior of \/3k\ 2 around this absolute minimum can easily be established. From a Taylor expansion, 
one obtains 



\Pk 



N 



(P-Pmin) 2 , |/3fe|" 



l(i 



J 2 iV 2 ($-$ min ) 2 



(31) 



For a crude order of magnitude estimate, one can develop the Bessel functions to first order in the small argument 
limit e Vm — 2\/2k c / K + <C 1 (more exactly, for de Sitter inflation and k c = Mp\, one has k c /K + ~ 0.06). This leads to 
J ~ 1/16+ 0{kl/K%) and N ~ ln 2 (\/2fc c /X + )/(47r 2 )(A' 2 /fc c 2 ) + 0(k°/Kl). Thus in order to avoid a back-reaction 
problem, the initial conditions p and $ must not differ too much from p m [ n and $ m i n which lead to (3k = (hence a 
zero amount of energy density created) . More precisely, the energy density produced is of the order of the background 
energy density, i.e. dfi w /dln(fc p h ys ) = 1, when p or $ respectively depart from the minimum by amounts Sp or <5$ 
given by: 



8p~ O 



H 



In" 



Pi 



pi 



H 



H 



In" 



pi 



Ah 



pi 



H 



(32) 



Here we assumed fc c = Mp\. Hence the corresponding fine-tuning of the initial conditions is of order H/Mp\ (if one 
assumes a uniform measure in p and $ in parameter space). 

One should note that the above constraint is valid for a given co-moving wave-number k and has been calculated at 
a time when k/a — k c . Since the fractional density parameter of quanta extracted out of the vacuum dr2 u /dln(fc p h ys ) 
must be smaller at all times during inflation, i.e. for a range of co-moving wave-numbers k since k and rj can be 
related by the above constraint k/a = k c , the above constraint on p m i n and $ m i n rather applies to a continuum of 
values of co-moving wave-numbers. In other words one does not have to fine-tune two parameters characterizing the 
initial data but a whole continuum of parameters, i.e. the functions Pmin(k) and $ m in(fc) themselves. The dependence 
in k of these functions is hidden in the argument of the Bessel and Neumann functions e am , since y m depends on k. 

B. Non conformal coupling: £ 7^ 1/6 

One can also perform a similar calculation of the Bogoliubov coefficient when the coupling is no longer conformal 
£ 7^ 1/6. In this case the calculation can be performed analytically for all background scale factors. For the sake of 
simplicity we choose minimal coupling £ = but this can be trivially expanded to various choices of the coupling to 
gravity, and does not modify the conclusions we derive below. 

If £ = 0, the evolution of the modes is dominated by a" /a in the tail, i.e. when w p i iys -C H (fc p h ys 3> ^c)j and the 
solution can be written as: 



(33) 



11 



where C+(k) and C_(fc) are two dimensionless k— dependent constants, and rji is some initial conformal time. One 
can check that this solution and the power-laws given in Eq. (0} are the same. Here one cannot choose the time 
of matching 7/ m to the WKB solution, since the matching has to be done when fi^ = 0, i.e. when ui 2 — a" /a or 
fcphys = K+- In the region w p h ys 3* H, i.e. for 77 3> rj mi the WKB approximation is valid as we have seen before, and 
the matching to the WKB form is justified at rj = rj m (k). For a given wave- number k, we are free to set 77; = r] m , since 
this amounts to a redefinition of the constants C+(fe) and C_(fc) by a function of k. The matching at rj m then gives: 



Oik 







wo 









Pk = 



y/AkuJk 



C + 7, + - - — 



C- 



(34) 



with the function 7^ = (u/ /2u>) + iui as above [see Eq. <|16|l ] and where all quantities in the above two equations are 
understood to be taken at 77 = r) m {k). In particular, at time ?7 m , uj 2 = a" /a = (1 + /3)a m iJ m //3. Since k p h ys l$> k c at 
rjra one can use the limiting form of <x> p h ys given in Eq. 10, hence: 



Ik 



a 
a 



a m H„ 




(35) 



The constants C + (fc) and CL(fe) are related to one another by the normalization of the mode functions: /-tfe/i£' 
MfeMfc = i, which gives: 



C+{k) = -p 



H m r sin $ ' 



(36) 



and as before we keep r = |C_(fc)| and <& = arg[C_(fc)] as the two independent parameters characterizing the choice 
of initial data. One finally derives the squared modulus of the Bogoliubov coefficient /3k as: 



\Pk\ 



P 2 



i 2 $ 4 



/ P K, 
P+lH n 



1 



4fc r 



1 + /3 
P 




cot$ 



(37) 



Since K + / H m is a large number, in the following we use the rescaled variable p = ry H / K + instead of r. Equation 11371) 
above is particularly attractive because it has exactly the same functional shape as Eq. (|18fl . It allows us to understand 
analytically the behavior of the amount of energy density produced in modes with fc p h ys ~ k c as a function of the 
initial data. The occupation number \(3k\ 2 has an absolute minimum located at: 




cos$ mill = 



1 + K+/(4fc c ) 



[l + K + /(4k c )} 2 + (l+p)/p 



(38) 



As before the occupation number vanishes exactly at this minimum, but the back-reaction problem cannot be avoided 
for generic initial conditions. In the present case it is not possible to make a sensible contour plot of d$l w /d ln(fc p h ys ) 
since this function changes by many orders of magnitude over very small intervals of p, $. Therefore, we take a 
conservative approach in which we calculate |/?fe(p)| 2 as a function of p for the values of <& that minimize this quantity 
at each p. We also evaluate |/3fc(<I>)| 2 as a function of $ for the values of p that minimize this quantity at each In 
other words, we solve c*$|/3/c| 2 = for $ as a function of p and d p \Pk\ 2 = for p as a function of $: 



$min(p) = tan 



f3 2 [l + K + /(4k c )] 2 + (l + f3)/f3 
p 2 l + K+/(Ak c ) 



4 

r mm 



sin 2 <f> 



4fc c 



1 + /3 



P 



(39) 



and plot dfL/dln(fc phys ) for \(3k[p,^min(p))\ 2 in Figs. and \Pk[Pmin(®), ®}\ 2 in Figs. El (for respectively de Sitter 
inflation and today). One clearly sees from these figures that dJ\,/d ln(/c p h ys ) ~ Mp^/H 2 for most values of which 
corresponds to our previous expectations, i.e. the amount of energy density created in quanta with k p hy S ~ k c ~ Mp\ 
is of order Mp V The behavior of \(3k\ 2 around the local minimum can be studied in the same way as before and one 
obtains: 



\h\ 2 {p) * ^ 



P 



\Pk\ 2 m 



1 + p 

P 
1 + P 



1/2 



(fl- 



ip ~ Pmin) , 



4fc r 



1 + /3 



(<b - <J> 



(* ~ *min)- 



(40) 



(41) 
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FIG. 3: Left panel: the solid line represents log 10 (|/3fc | 2 /|/3™ ax | 2 ) = log 10 [dO„/d ln(fcp hys )] plotted as a function of the rescaled 
variable p characterizing the initial data. This plot corresponds to an inflationary era with a de Sitter metric, and Hubble 
parameter ~ 10 _6 M P i. The other parameter of initial data is $ = <E> m m- Allowed regions correspond to log 10 (|/3fc| 2 /|/9™ ax | 2 ) < 

0, and are peaked around a particular value of p. The minimum is in fact \/3k\ 2 = 0, corresponding to log 10 (|/3fc| 2 /|/3fe lax | 2 ) = 
— oo, but cannot be seen in the figure due to insufficient resolution. The dotted line provides a continuation of the numerical 
result to the analytical value at that point. In most of parameter space, the energy density is too large by ~10 orders of 
magnitude. Right panel: Same as left panel for Ho ~ 10 _61 Mpi in a matter dominated era. In nearly all of parameter space 
the energy density is too large by ~ 122 orders of magnitude . 





FIG. 4: Left panel: log 10 (|/3fc| 2 /|/3™ ax | 2 ) = [dfi„/d ln(fe p hy S )] plotted as a function of the phase characterizing the initial data. 
This plot corresponds to an inflationary era with a de Sitter metric, and Hubble parameter ~ 10 _6 Mpi. As indicated in the 
text, the other parameter of initial data (a modulus) has been eliminated for the phase $ by minimizing \/3 k \ 2 . Allowed regions 
correspond to log 10 (|/3fc| 2 /|/3™ ax | 2 ) < 0, and are seen to be peaked around a particular value of The plot is symmetric in 
the interchange <E> <-» $ + 7r. Note that the ordinate scale is in log 10 : in most of parameter space, the energy density is too large 
by ~10 orders of magnitude. Right panel: Same as left panel for H ~ 10 _61 M P i in a matter dominated era. In nearly all of 
parameter space the energy density is too large by ~ 122 orders of magnitude, i.e. dfi^/d m(fc p hy S ) ~ Afpi- 



If we write the value |/3fc|max such that df2 w /dln(fc p hy S ) = 1 (i.e. similar amount of energy created in quanta with 
fcphys ~ k c than in the background), then this value is reached if p and <I> departs from p mm and ^min by an amount 
Sp, 5$, with: 



5p~ O 



H 



pi 



5$ - o 



H 

mT 



•In" 



Pi 



M P1 
H 



(42) 



In other words, the fine-tuning necessary to avoid back- reaction is of order ~ <5$ ~ Sp ~ H/Mp\ (neglecting the In 
term). During inflation H/Mp\ ~ 10~ 6 , and today H/Mp\ ~ 10 -61 . Although the fine-tuning can be considered as 
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not too severe during inflation, we see that, during the matter epoch, it is as severe as the usual fine-tuning of the 
cosmological constant problem. Therefore, the scenario of Refs. 0,0,13 does not improve the situation. 

At thispoint, it should be emphasized that the above problem is a generic feature of the dispersion relation used 
in Refs. U H 0. In order for the trans-Planckian effects to modify the power spectrum of the fluctuations, the 
physical modes of interest must spend some time in a region where the WKB approximation is violated. As already 
mentioned, this implies production of particles and the energy density associated to these particles must not exceed 
the background energy density. This implies some constraints on the occupation number \Pk\ 2 - It has been shown 
in Ref. |8| that these constraints are quite stringent if the production is taking place today. Usually, these tight 
constraints can be avoided by requiring that the violation only occurs during inflation where the problem is less 
severe. This can be achieved if the dispersion relation is such that w p h ys 3* H for all trans-Planckian wave-numbers, 
with H the Hubble constant some unspecified time after inflation. Such an example is provided in Fig. 2 of Ref. [||. 
Then as was argued in the discussion following Eq. 10 the WKB approximation should be valid at all times after 
inflation and adiabaticity restored. Here however, since w p hy S — > as fc p h ys — * +00, there is at all times a region where 
the WKB is violated. Therefore the class of dispersion relations used in Refs. 0, Li L| suffer in a generic way from 
the problem discussed in Ref. 0- 

Furthermore, we also note as previously that the above calculation of the fine-tuning holds for a given co-moving 
wave-number k. Similar constraints apply for other wave- numbers but the values of $ m i n and p m { n are shifted from 
the above [notably because of the choice rji = rj m (k)]. Therefore one must not only pick one right value of $ to one 
part in ~ Mp\/H, but a whole continuum of values of <5(fc) such that back-reaction is avoided for all of these values. 



IV. EQUATION OF STATE 

Up to now we have argued that: (i) in the scenario proposed in Refs. P-LllQ there is no preferred initial vacuum 
state, hence all conclusions drawn depend on the ad- hoc choice of initial data; (it) for arbitrary values of the two 
parameters characterizing this choice of initial data one finds that energy in excess of the background energy density 
is produced due to the non adiabatic evolution of modes in the tail. 

In a separate publication, we have constructed an effective energy-momentum tensor for scalar field theories with 
non-linear dispersion relations, and we have shown that dispersion relations of the form of that proposed in Refs. 0, 
0, 13 generically led to the wrong equation of state. This finding has been challenged by Bastero-Gil and Mersini 
recently, who argue that the energy-momentum tensor we have constructed is ill-defined as it (supposedly) is not 
divergenceless. 

Explicitly, in Ref. [g , the vev for the energy density and pressure are given by 



(p) 
(P) 



1 



4ir 2 a 4 
1 

47r 2 a 4 



dkk 2 



dkk 2 



(?)' 



■to 2 (k) \p k \ 



2 2 dw 2 

3 dfc^ 



u 2 ) |/ife| 2 



(43) 
(44) 



and the integrals extend from to +00. Reference Q claims that (p)' + 3H(p + p) 7^ 0, with "H 
noting that the co-moving frequency w(k) = aa> p h ys (&phys), a straightforward calculation gives: 



a' la. However, 



(pY 



Air 2 a 



dkk' 



4 - Hp! k + { H 2 - — ) /i fe 



/4 - Kp* k + [ n - — 1 p* k 



04' - Kpi) + (p' k - Hp k ) 

-AH - Hp k f + oj 2 (k) (p' k p* k + p k p* k ') - 2Hk 2 ^\p k \ 2 - 2Hu 2 (k)\p k \ 2 } 

4^4 J dfcfc2 | (Vfc - i^k - U P*k) + {Pk - HPk) [pf - ~P*k) - 67i \p'k - KpkX 

du 2 

+u 2 (k) [fi k {p* k ' - Hpl) + pi (p' k - Hp k )] - mk 2 — \p k 



^ f dkk 2 (-6H\p' k -Hp k f -2Hk 2 ^\p k \ 2 



(45) 
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where the field equation p'l — (a" /a)pk = — w has been used in the last step. Since 

ZH{P + P) = / dkk2 (enlfii-Hfikf + 2Hk 2 ^\p k \ 2 ^ , (46) 

the energy conservation condition (p)' + 3H(p + p) = is trivially satisfied. 

We take advantage of this Section to point out that the energy-momentum tensor we have constructed in a previous 
publication 6] is well behaved and the construction is entirely consistent. In effect we constructed a generally covariant 
Lagrangian for a scalar field with a non-linear dispersion relation. The breaking of the Lorentz invariance is explicited 
by introducing a dynamical four-vector in the Lagrangian whose role is to define the preferred rest frame while 
preserving general covariance at the same time, following previous work by Jacobson and Mattingly jl3l | (see also 
references therein). The energy-momentum tensor derived by varying the action with respect to the metric is the sum 
of the energy- momentum tensor of the scalar field and that of it' 1 . We have restricted ourselves to FLRW space-times 
by choosing u M = (—1,0,0,0). This approach is consistent as satisfies its field equation [see Eq. (15) in Ref. 0], 
and the scalar field also satisfies its field equation. In FRLW space-time, the energy-momentum tensor of the four- 
vector (i.e. the part which depends only on u^) vanishes since is constant 0, hence the remainder is the scalar 
field energy-momentum tensor, which is conserved as shown above. We consider this scalar field as a test field, which 
we quantize on the curved background. The dynamics of the background can be specified by adding matter fields 
to the action without modifying our approach. Indeed the matter field energy-momentum tensor will be separately 
conserved. Therefore our earlier criticisms on the equation of state of the trans-Planckian modes apply. 

We also note that a recent paper by Frampton [T(j argues that by adding higher order terms of the form u^V^Un 
in the effective Lagrangian, one can obtain a correct equation of state, i.e. similar to that of vacuum. Here T> 2 denotes 
a three-dimensional Laplacian expressed in a generally covariant way on space-like hyper-surfaces (see Refs. 0,113 for 
explicit definitions). The choice — (— 1, 0, 0, 0) then reduces this term to the usual three-dimensional Laplacian on 
constant time hyper-surfaces for FLRW space-times. However, as should be obvious, the only terms that can enter the 
energy-momentum tensor via this new term always contain derivatives of the form T> 2n u il when n > 2. Those terms 
vanish when one picks to be constant as above. For n = 1, i.e. the lowest order term of the expansion, one can 
check that the variation of 2? 2 it M with respect to the metric or its first derivative always induces a term proportional 
to a derivative of u^ 1 that is either spatial or temporal. This calculation can be found in Eqs. (B3) and (B4) of Ref. @ 
(the substitution <j> — > in this calculation can be made without changing the result since the equation is written in 
a generally covariant way). Therefore we conclude that the addition of terms proposed by Frampton cannot lead to 
any difference with respect to the energy-momentum tensor we proposed earlier. Hence such terms cannot account 
straightforwardly for a vacuum-like equation of state in this scenario. 

The difference with the conclusion of Ref. ^(j probably lies in the choice of normalization of u M : Ref. [lj} claims 
that, in a cosmological context, u^u^ — a" 2 = (a is the scale factor). However, u^u^ is a scalar, while a~ 2 is the 
component of a rank two tensor so that this choice of normalization is manifestly not covariant. Of course, it can be 
made covariant at the expense of the introduction of a second vector field but this does not seem to be the case in 
Ref. This implies that the Lagrangian proposed by Frampton is generally not covariant. Of course the correct 
normalization for in a cosmological context is always u^u^ — 1 = (up to the sign convention), and u^ L can be 
written as (— 1/a, 0,0,0) if the FRLW metric is written in terms of conformal time, or (—1,0,0,0) if the metric is 
written in terms of cosmic time. 

To finish, let us also note that it is claimed in Ref. pj that "the tail modes are still frozen at present... Thus 
the energy of the tail is a contribution to the dark energy of the Universe: up to the present it has the equation of 
state of a cosmological constant term" . Here, that the tail modes are frozen means (pk/a)' ~ 0, the solution to the 
field equation when the term a" /a dominates (for £ = 0). But there is no logical relationship between these modes 
being frozen and the equation of state being that of a cosmological constant. As a matter of fact the equation of 
state of the modes of wavelengths larger than the horizon size for a scalar field with a linear dispersion relation takes 
the form p = —p/3. The derivation of the equation of state requires to define correctly a stress-energy tensor when 
such non-linear dispersion relations are taken into account (and thus when Lorentz invariance is broken) as we did 
in Ref. 6] . Note also that there is a contradiction between the above claim of Ref. pj and the whole content of the 
recent paper by Bastero-Gil and Mersini |3| which calculates the equation of state of the trans-Planckian modes and 
finds that it is not the equation of state of the vacuum but that it approaches it only at late times. 



V. CONCLUSIONS 

In this article, we have examined in detail the scenario proposed in Refs. QHHQ which attributes the dark energy 
to the properties of a scalar field with a dispersion relation that decreases exponentially with trans-Planckian wave- 
number fcphys k c ~ Mpi. We have demonstrated that this mechanism does not work, mainly for two reasons. 
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(i) The mode function of the scalar field does not behave as a plane wave as r\ — > — oo in the so-called "tail" (i.e., 
the part of the non-linear dispersion relation where Wph ys <C H and fc p h ys 3> Mpi). This implies that there is 
no definite prescription for constructing a well-defined initial vacuum state, hence that the choice of initial data is 
entirely arbitrary. This situation is similar to the problem of setting initial data for cosmological perturbations in the 
absence of an accelerated expansion era (but with linear dispersion relations), in which case the data would have to 
be specified on super-horizon scales where the mode function does not oscillate and is frozen by the expansion. In 
the scenario of Refs. Q, 0, the WKB approximation is not valid at all times for modes in the tail. This explains 
that the notion of adiabatic vacuum cannot be used to set up initial data and the initial state proposed in Ref. 
is thus ad-hoc. It also brings us to the second objection against this scenario: (ii) since all modes originate in the 
tail of the dispersion relation, the breakdown of the WKB approximation for a given physical wave-number at early 
times implies the continuous production of a substantial number of quanta with physical wave- numbers > k c . The 
breakdown of the WKB approximation can indeed be seen as the signal of a strongly non-adiabatic evolution which 
is generically associated with particle production in expanding space-times. We have calculated the amount of energy 
density produced in modes of physical wave-number k c as a function of the two parameters that characterize the 
(arbitrary) choice of initial data. We have shown that this energy density is generically of order Mp V The production 
of energy density in quanta with wave-numbers ~ Mp\ in excess of the background energy density ~ H implies 
the breakdown of the perturbative semi-classical framework used for the calculation, and renders all claims irrelevant. 
There is a small region of parameter space in which this energy density can be tuned down to zero, but the fine-tuning 
in the choice of initial data is of order H/Mpi. Such a fine-tuning at the time of inflation is of order ~ 10~ 6 and 
is probably acceptable, but today, it is of the same order as the celebrated fine-tuning of the cosmological constant 
problem. 
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